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ABSTRACT 

We have made a Monte Carlo simulation of the intergalactic absorption in 
order to model the Lyman continuum absorption, which is required to esti- 
Q^ ■ mate the escape fraction of the Lyman continuum from distant galaxies. To 

(N 
^' 

' intergalactic absorbers which reproduces recent observational statistics of the 

oo : . . , , 

, Lyman a forest, Lyman limit systems (LLSs), and damped Lyman a systems 

(DLAs) simultaneously. In particular, we assume a common functional form 
' of the number evolution along the redshift for all types of absorbers. The 

' Lyman series transmissions in our simulation reproduce the observed redshift 

evolution of the transmissions excellently, and the Lyman continuum trans- 
mission also agrees with an observed estimation which is still quite rare in the 
literature. The probability distribution of the Lyman a opacity in our simu- 
lation is log-normal with a tail towards a large opacity. This tail is produced 
by DLAs. The probability distribution of the Lyman continuum opacity in 
our simulation also show a broad tail towards a large opacity. This tail is 
produced by LLSs. Because of the rarity of LLSs, we have a chance to have a 
clean line of sight in the Lyman continuum even for z ^ A with a probability 
of about 20%. Our simulation expects a good correlation between the Lyman 
continuum opacity and the Lyman a opacity, which may be useful to estimate 
the former from the latter for an individual line of sight. 
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1 INTRODUCTION 

Neutral hydrogen rem aining in the intergalact ic medium (IGM) absorbs the radiation from 
distant sources (e.g., Gunn fc Peterson IQGsl ). This intergalactic absorption seems to be 



caused by numerous discrete systems, which are called absorbers in this paper, on the line of 
sight. The absorbers probably consist of various types of objects: the outer edge of galaxies, 
halo gas, diffuse medium in the interg alactic space, etc. They probably connect with each 



other and form the "cosmic web" (e.g. 



Ranch 



19981). 



Because of such clumpy nature of the IGM, the intergalactic absorption fluctuates among 
the lines of sight. Thus, we cannot predict the amount of the absorption for a certain 



dard deviation for many lines o 


" sight 


in a statistical 


Giallongo & Cristiani 


1990; 


Zuo 


1993; 


Madau 


1995 


)• 



sense (e.g.. 


M0h 


er & Jakobsen 


1990; 


Bershadv et al. 


( 


1999 


) showed that a 



Monte Carlo simulation is very powerful to discuss not only the mean absorption but also the 
dispersion correctly. Here we present such a simulation. Because many observational results 
of the IGM absorbers have been pubhshed since 1999, we update the absorbers' statistics 
to be input ii ito the s ir aulati on. 

Recently, iMeiksinI (120061 ) presented a new mean curve of the intergalactic absorption 



based on the density distributio n produced by a cosmological s imulation. However, he dis- 



cussed mean transmission only. 



Tepper-Garcia &: Fritzd (120081 ) have done a Monte Carlo 



simulation and discuss the distribution of the absorption amount. However, they concen- 
trate on only Lyman a absorption. In this paper, we will examine the stochastic nature not 
only of the Lyman series absorption but also of the Lyman continuum one. 

This is motivated by the recent observational attempts for determining the escape 



fraction of the Lyman con 


;inuum, especially the hydrogen ionising photons, from c 


istant 


galaxies (e.g.. 


Steidel et al. 


2001 


Gial] 


oneo et al. '. 


2OO2I; 


Vlalkan. Webb. & KonoDackev 


2003; 


Inoue et al. 


2005; 


Shaplev et al. 


2006 


Siana et al 


. 200? 


). In the process towards the escape 



fraction, we need a correction of the intergalactic absorption for the Lyman continuum of the 
distant galaxies. Since the Lyman continuum absorption by the IGM has not been measured 
well observationally, a model predicting the absorption amount of the Lyman continuum is 
required. 
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The effective optical depth throu gh a clumpy IGM for the frequency in the rest-frame 



of a source at the redshift zs is (e.g. 



Paresce et al. 



1980) 



dz 







Ni ozoNui 



(1) 



where d'^N'/dzdNm is the number of absorbers along the line of sight per unit redshift z 
interval and per unit HI column density A^hi interval, and = crHi(i^s(l + ^)/(l + zs))Nm 
is the optical depth of an absorber with A^hi at z, with being the HI cross section o"hi(z^) at 
the frequency u in the absorber's rest-frame. If we assume the column density distribution 
of the absorbers as a power -law with an index —(3 independent of the redshift, which is well 



found in observations (e.g.. iTytler 



19871 ). we have 



HI 



e ■^"M OC 



2-/3 
HI 

1-/3 
HI 



(TcI < 1) 
(TcI > i: 



(2) 



Since observations show /5 ^ 1.5 (e.g. 
made by absorbers with Td ~ 1 (e.g.. 



Tvtlei 



19871). the maximum contributi on to is 



Mciller fc Jakobsen 



1990 



Meiksin 



20061). Thus, the 



absorption of the Lyman continuum is mainly caused by the Lyman limit systems (LLSs) 
with A^Hi ~ 10^'^ cm~^, not by the Lyman a forest (LAF) with A^hi ~ 10^^ cm^'^. Therefore, 
we take a great care of the treatment of LLSs in order to properly compute the Lyman 
continuum absorption by the IGM. Moreover, the Lyman continuum absorption should 
be very stochastic because LLSs are relatively rare. Therefore, we also need to model the 
dispersion of the absorption correctly. 

Previous works about the intergalactic absorpt i on as sumed different number evolutions 
along the redshift for the LAF and LLSs. iMadaul (119951 ) assumed dj\ f /dz oc (1 -|- zY/'^^ fo r 



the LAF (Murdoch et al 



Bershady et al 



19861) and dAf/dz oc (1 + z)^-^^ for LLSs dSargent et al 



19891). 



(119991) ass u med a different set of the LAF number evolut io n and th e same 



LLS evolution as 



Madau (1995) (see also Tepper-Garcfa fc Fritze 20081 ). MeiksiJ (2006) 



assumed the number evolution extracted from a cosraologic al simulatio n for the LA 
dAf/dz oc (1 + zY'^ for LLSs (IStengler-Larrea et al.lll995l ). However 
showed in their appendix that recent LLSs observations reported by 



Inoue et al. 



Peroux et al 



and 



(2005) 



(l2003h 



can be reproduced by a common number evolution function for both types of absorbers. We 
take this common function scenario in this paper because of its simplicity. 

The rest of this paper consists of 5 sections; in section 2, we present an empirical dis- 
tribution function of the intergalactic absorbers and show that the function reproduces all 
the observational data simultaneously. In section 3, we describe the procedure of our Monte 
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Carlo simulation. In section 4, we describe the method for calculating the IGM transmission, 
compare our transmission model with others, show that our model reproduces the observed 
Lyman series transmissions very well, predict the Lyman continuum transmission, discuss 
the distribution functions of the intergalactic opacities of the Lyman continuum as well as 
Lyman a line, and present a method for estimating the Lyman continuum opacity from the 
Lyman a one. In section 5, we discuss the detectability of the Lyman continuum from distant 
galaxies with our Monte Carlo results. The final section is devoted to our conclusions. 



2 DISTRIBUTION FUNCTION OF INTERGALACTIC ABSORBERS 



In this section, we present an empirical distribution function of the intergalactic absorbers 
which is input into the Monte Carlo simulation in the later sections. We consider a func- 
tional form with several parameters which are determined so as to reproduce the observed 
redshift, column density, and Doppler parameter distributions of the intergalactic absorbers. 
Especially, we assume a common number evolution function along the redshift for all types 
of absorbers. 

An absorber have three physical quantities characterising itself: redshift z, HI column 
density A^hi? and Doppler parameter b. Here, we assume that these three quantities are 
independent of each other. This is just for simplicity although a slight correlation between 



Kim et al. 



200ll ). Then, we express the number of absorbers per 



A^Hi and b has found (e.g. 
unit volume of the three parameters' space as 

= (3) 

Note that we have assumed a common functional form of the number evolution along the 
redshift for all types of absorbers as f{z). In the following, we determine the functional forms 
of f{z), g{Nm), and h{b) so as to reproduce the observed distribution functions. Since the 
observed statistics of absorbers are not very tight yet, we do it by visual inspection instead 
of doing it by a rigorous way like a likelihood method. 

Figure 1 shows the column density distribution. The vertical axis is the number of ab- 
sorbers per unit column density and per uni t "absorption distance" which w as introduced 
to remove the effect of the Hubble expansion (iBahcall fc Peebleslll969l ). While iTytlerl ( 119871 ) 



proposed a single power-law for the column density distribution, t 
double power-law with a break at ~ 10^'' cm^^ as found in Figure 1. 



le recent data seems a 



Prochaska et al 



fl2005h 



also suggest a similar break around 10 cm . For the highest column density, however. 
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logio ^Hi [cm"^] 

Figure 1. Number of the intergalactic absorbers per unit column density per unit "absorption distance" along an average line 
of sight as a function of the column density. The shaded area is the observed range of the Lyman a forest (LAF) at z^bs = 0.5- 
1.9 by Janknecht et al. (2006). The filled circles are the data of the LAF at z^bs = 1.5-4.0 based on Kim et al. (2002). The 
diamonds, triangles, and squares are the observed data of sub-damped Lyman a systems (sub-DLAs) and damped Lyman a 
systems (DLAs) by Peroux et al. (2005; Zabs = 1.8-5.0), O'Meara et al. (2007; ^abs = 1.8-4.2), and Prochaska et al. (2005; 
2^abs = 2.2—5.5), respectively. All the data are scaled for a flat A cosmology with S7a = 0.7. The histograms are the column 
density distributions of absorbers generated by our Monte Carlo simulation: the absorber's redshift ^^bs = 0-2 (bottom), 2-4 
(middle), and 4-6 (top). The vertical shifts of the histograms are the result of the number density evolution along the redshift. 
The two solid lines are for the reference of the power-law index. 




0.01 U I > > I 
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Figure 2. Number of the intergalactic absorbers per unit redshift along an average line of sight as a function of the absorbers' 
redshift. The shaded area is the observed range for absorbers with log]^g(A^Hl/cm~^) > 13.6 (LAF) taken from Weymann 
et al. (1998), Kim et al. (2001), and Janknecht et al. (2006). The filled circles, triangles, and squares are the observed data 
of absorbers with log]^o(^Hl/cm~^) > 17.2 (LLS) taken from Peroux et al. (2005), > 19.0 (sub-DLA) taken from Peroux 
et al. (2005), and > 20.3 (DLA) taken from Rao et al. (2006), respectively. The number density evolutions of absorbers 
generated by our Monte Carlo simulation are shown by different symbols depending on the column density range: crosses 
(logio(7VHi/cm~^) > 13.6), open circles (logjQ(7VHi/cm~^) > 17.2), open triangles (logj^o(^Hl/cm~^) > 19.0), and open 
squares (logjg(7VHi/cm~^) > 20.3). They excellently trace the assumed distribution functions of absorbers given in equations 
(3)-(6), but we see statistical fluctuations in large column density cases due to the small number of such absorbers. 
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there may be a more rapid decline of t he number of absorbers, suggesting; a Sche chter hke 



function for this column density range (iPeroux et al. 



2003 



Prochaska et al. 



20051 ). Here we 



adopt a double power-law for simphcity as 

-A 

-132 



HI 



B 



(4) 



(iVc ^ iVffl ^ iV,) 

with /3i = 1.6 and (32 = 1.3 obtained by eye. The break column density Nc is assumed to 
be 1.6 X 10^^ cm^^ which is the usual criterion for LLSs, i.e. the absorber with a column 
density larger than this value is optically thick for the Lyman limit photon. A physical 
explanation for this break, which gives more number of higher column density absorbers, 
is a self-shielding effect; the neutral fraction inside an optically thick absorber (i.e. LLSs) 
could be kept higher. The normalisation B is determined by / giN-aijdNin = 1. 

Figure 2 shows the number evolution along the redshift of the LAF, LLSs , sub-DLAs, and 



DLA s. The LAF number evolution (the shaded area) shows a break at 2; ~ 1 (jWeymann et al. 



19981). 



Dave et al. 



(1l999l ) have explained that this break is due to a steep decline of the 
ionising background intensity from z ~ 1 to 0. Thus, we assu me a break at 2; ~ 1 in the 
functional form of f{z). On the other hand, iFan et al.l (120061 ) suggest a steepening of the 
Lyman a opacity ev olution at z_ >_4, which means more absorbers than expected at high 

(l2006l ). thus, we assume another break at z ~ 4 in /(z). 



redshift. As done in 
Therefore, we assume 

' l+z 

f{z)=A< 



Inoue et al. 



71 



Vl+zi 
L \l+zij Vl+22/ 



(5) 



{0< z ^ zi) 

{Zi < Z ^ Z2) . 

{Z2 < Z) 

We adopt zi = 1.2, 71 = 0.2, and 72 = 2.5 which are determined by eye to fit the observed 
number evolutions along the redshift of the LAF in Figure 2. We also adopt the second break 
in f{z) a t Z2 = 4.0 with 73 = 4.0 to reproduce a rapid increase of the Lyman a opacity 
found by iFan et al.l (120061 ). With these parameters, the observed number evolutions of all 
types of absorbers are reproduced simultaneously. Note that we have assumed g^Nm) to be 
a double power-law as equation (4). If we assumed 5'(A^hi) to be a single power-law, we could 
not reproduce the number evolution functions of all types of absorbers simultaneously. The 
normalisation A is the total number of the IGM absorbers at z = Zi with a column density 
Ni ^ A^Hi ^ A^u- We adopt A = 400 with Ni = 10^^ cm^^ and = 10^^ cm'^ to match 
with the observed number of the LAF in Figure 2. 
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Table 1. Parameters for the distribution function of intergalactic absorbers. 



Parameters 




Adopted value 


71 




0.2 


72 




2.5 


73 




4.0 


zi 




1.2 


Z2 




4.0 


A 




400 


logioW/cm" 




12.0 


logio(^u/cm 




22.0 


logio(^c/cm- 




17.2 






1.6 


& 




1.3 


fecr/km s~^ 




23.0 



For the Doppler parameter distribu tion function h{b ) , we assume the functional form 



with a single parameter b^ suggested by 



Hui fc Rutledgj fllQQQf ): 



hib) 



55 



(6) 



Janknecht et al 



(120061). We note 



We adopt b^ = 23 km s based on the measurements by 

POO 

that h(b) in equation (6) is already normalised as / h(b)db = 1. 

Jo 

In Figures 1 and 2, we show comparisons of the empirical distribution functions presented 
in equations (4) and (5) with the observed functions. For our distribution functions, we 
show the results generated by our Monte Carlo simulation, whose procedures are described 
in the next section. We have confirmed that our Monte Carlo simulation reproduces the 
input distribution functions excellently. Figure 1 shows that our column density distribution 
g{N-iii) is very consistent with the observations. Note that the vertical axis in the panel is not 
exactly same as g{Ni{i) but the column density distribution averaged over a redshift range 
(and divided by the "absorption distance"). Figure 2 shows that our redshift evolution f{z) 
nicely reproduces the observed redshift evolution for all types of absorbers simultaneously. 
Note that the vertical axis in the panel is not exactly same as f{z) but the number density 
of absorbers with column densities larger than a certain value. Although the disagreement 
of DLAs at 2 ~ 1 (but less than a factor of two) might suggest a redshift dependence of 
g{Nm), we avoid this complexity. The adopted values for the parameters are summarised in 
Table 1. 

Figure 3 is the close-up of the number evolution of LLSs along the redshift. As described 
in section 1, we need to reproduce the observed number of LLSs for a rigid prediction of 
the Lyman continuum absorpti on by the IGM. Figu re 3 ensures a good agreement between 



our model and observations by 



Peroux et al. 



20051 ). On the other hand, the regression line 
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Figure 3. Number of LLSs per unit redshift along an average line of sight as a function of the LLSs' redshift. The filled circles 
are the recent observed data by Peroux et al. (2005). The open circles are our model. The dashed line is the regression line by 
Stengler-Larrea et al. (1995), which is assumed in Meiksin (2006). The dotted line is the regression line by Sargent et al. (1989), 
which is assumed in Madau (1995), Bershady et al. (1999), and Tepper-Garcia & Fritze (2008). 



sussest e d by 
MadaJ ( 19951 ) 



Sargent et al. 



Bershady et al. 



198^ sho wn as the dotted line in Figure 3, wh ich is assumed in 



(119991 ) , and iTepper-Garcfa fc Fritzd (120081 ) , does not agree 



well with the recent data. T he dashed line in Fi gure 3 is the r egression line suggested by 



Stengler-Larrea et al 



(119951 ) which is assumed in iMeiksinI (120061 ) . This case agrees with the 
recent data. However, in this case, we would need different functional forms of f{z) for the 
LAF and for LLSs to fit all the data shown in Figure 2 simultaneously. 



3 MONTE CARLO PROCEDURE 

Based on the assumed distribution function described in the previous section, we gen- 
erate a large number of absorbers along lines of sight by a Monte Carlo method. The 
random number generator use d in this paper is the "Mersenne Twister" developed by 
Matsumoto fc Nishimural ( 19981 ) which is a very fast random number generator with a very 
good statistical property, l 

We assume the encounter of an absorber on a line of sight to be a Poisson process. In 



other words, we neglect the effect o : 



found (lOstriker. Bajtlik. fc Duncan 



clustering of absorbers although it has been already 



19881 ). We should examine the effect in future. For a 



Poisson process, if we have an absorber at z, the probability encountering the next absorber 
at z + Az is 



^ The "Mersenne Twister" is proved to have the period of 2^^®^'' — 1 and the 623-dimension equidistribution property. Its 
C-code is public at |http://www. math. sci.hiroshima-u. ac.jp/~m-mat/MT/emt. html| 
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piAz;z) = fiz)e~f(^^'^\ (7) 

because the mean redshift interval of two absorbers at z is just the reciprocal of the redshift 
distribution function f{z) of equation (5). 

For a line of sight, we start from z = and determine the redshift of the first absorber 
by drawing a random number based on equation (7), and we determine the column density 
and the Doppler parameter of this absorber by drawing other two random numbers based on 
g{Nm) and h{b) of equations (4) and (6). We note that 5'(A^hi) and h{b) themselves are the 
probability distribution functions. Then, we determine the redshift, the column density, and 
the Doppler parameter of the next absorber likewise. These procedures are repeated until 
the absorber's redshift exceeds z = 6. Typically, we generate about 18,000 absorbers for one 
line of sight which depends on the lower limit of the column density A^i. For a smaller A^i, we 
need much more absorbers which is time-consuming. As shown in equation (2), absorbers 
with a column density less than A^i = 10^^ cm~^ do not contribute to the opacity very much 
because their optical depths are much less than unity. We generate 10,000 lines of sight for 
enough statistics. 



4 INTERGALACTIC TRANSMISSIONS 

Based on the intergalactic absorbers along lines of sight generated by the Monte Carlo 
procedures described in the previous section, we calculate transmissions along the lines of 
sight. 

Suppose an absorber with the redshift z, the column density Nm, and the Doppler 
parameter b. The optical depth of this absorber for the frequency u in the absorber's rest- 
frame is 

T{u)=NuJai^c{^) + y2a,{u)] , (8) 



where (Tlc is the cross section for the Lyman continuum, and cTj is the cross section of the 
i-th Lyman series line. The cross section for the Lyman continuum is approximated to 



where the cross section at the Lyman limit z/ll is ctll = 6.30 x 10 cm^ ( 10sterbrocklll989l ) . 
When u < ui^i^, c'"lc(^) = 0. The cross sections for the Lyman series lines are approximated 
to 
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(10) 



meCZ/D 

where mg is the electron mass, e is the electron charge, c is the speed of light, fi is the oscilla- 
tor strength of the z-th line, z/^ = Vi{h/c) is the Doppler width with the line centre frequency 
Vj, and ct)i{u) is the lin e profile function which is calculated by an analytic approximation by 



Tepper-Garcfal (120061 ). The hydrogen atomic d ata of Pj, fj, and th e damping constant Fj for 



Wiese et al. 



(11966h . 



the line profile up to 40-th line are taken from 

Since we consider only Hi absorbers in this paper, we calculate the wavelength range 
not affected by helium: 700-1300 A in the rest-frame of the source redshift. Note that we 
should consider wavelength enough longer than that of the Lyman a line (1216 A) in order 
to treat a broad wing of damped Lyman a lines. We should also calculate transmissions with 
an enough high wavelength resolution. We have found experimentally that typical errors in 
transmissions within a wavelength range relative to those calculated with 0.01 A resolution 
in the source rest-frame are 10% for 1 A resolution and 1% for 0.1 A resolution^ Thus, we 
have calculated all the results in this paper with 0.1 A resolution in the source rest-frame. 



4.1 Example and average transmissions 

Figure 4 shows the IGM transmissions along an example line of sight. We find a lot of 
narrow absorption lines as the LAF. In the Lyman continuum region, we find several step- 
wise depressions caused by LLSs. For example, in the panel (b) (^s = 3 case), we find 
significant transmissions below the Lyman limit up to 864 A at which a sharp depression 
appears because of a LLS at z = 2.78 which has A^hi = 6.5 x 10^'' cm^^. This means that 
the Lyman continuum absorption by the IGM is very stochastic because it is controlled by 
the presence of relatively rare LLSs. This point is illustrated in Figure 5. If we do not have 
a LLS (or DLA) near the source, we can expect a significant transmission even far below 
the source Lyman limit (panel [a]). However, If we have such a LLS, the transmission is 
suddenly cut down at the corresponding wavelength (panel [b]). 

On the other hand, such stochastic behaviour disappears in average transmissions shown 
in Figure 6, where we also compare our models with those by other authors. The solid 
line in each panel is our Monte Carlo model. We have averaged 10,000 lines of sight in each 

^ We should note that relative errors are larger for lower transmissions. For a very small transmission, say < 10""^, the relative 
error becomes about 10% or more in some cases even with 0.1 A resolution. Such cases appear in the Lyman continuum 
transmissions for the source redshift 23 > 5. However, the fraction of the cases is ~ 10%. 
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700 800 900 1000 1100 1200 1300 

Rest wavelength [A] 

Figure 4. Examples of the IGM transmissions along a line of sight. The horizontal axis is the wavelength in the source 
rest-frame. The source redshifts are noted in the panels. 
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Figure 5. Example IGM transmissions of the source redshift zq = 3 for different lines of sight. 
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Figure 6. Average IGM transmissions. The horizontal axis is the wavelength in the source rest-frame. The source redshifts 
are noted in the panels. The solid lines are our Monte Carlo results; 10,000 lines of sight are averaged in each panel. The 
dot-dashed lines are the mean transmission models of Madau (1995). The dotted lines are the mean transmission models of 
Meiksin (2006) but the updated version. The dashed line in the panel (c) is the average transmission of the Monte Carlo 
simulation by Bershady et al. (1999) (MC-Kim model). 

panel. Thus, statistical variations are su ppressed in a very sma ll level. The dashed line in the 



panel (b) is the Monte Carlo model by 
the adopted distrib 
agreement between 



Bershady et al. 



(119991 ) (MC-Kim model). Although 



the adopted distrib ution function o f the absorbers is different from ours, we find a good 

(119991 ) and us for the source redshift zs = 3. This is 



Bershady et al. 



caused by a coincidence of the number density of the absorbers around 2; = 3 in our model 
with that in their model. For example, the coincidence of the number of LLSs around z = 3 
in our model with that in their model shown in Figure 3 makes a very good agreement of 
both transmissions below the Lyman limit. However, we may expect that some differences 
appear for other redshifts because of differences of the number density of the absorbers in 
other redshifts. 

As the dotted lines in Figure 6, we show the mean transmission models by 



Meiksin 



(120061 ). This is not the original version but the updated one where the treatment of the 
Lyman continuum absorption by the LAF was revi s ed (A . Meiksin, private communica- 
tion). The agreement between the updated iMeiksinI (120061 ) and us is excellent for all the 



redshift in the Lyman series regime, 
smaller transmissions than those of 
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n the Lyman c ontinuum regime, our simulation shows 



Meiksinl (120061) . This i s beca use the number of LLSs 



in our model is almost a. 



by IStengler-Larrea et al. 



Meiksin 



LLSs in 
is assumed in 
than ours and 



ways larger than that in 



1995 



Meiksinl ([20061) (i.e. the regression line 



see Figure 3). For example, ar ound z = 3, t h e num ber of 



Bershadv et al. 



Bershady et al. 



Sargent et a. 



19891 ) which 



20061 ) is abou t two -thirds of ours and that of 

(11999 ). so that the transmission of iMeiksinI ( l2006l ) is larger 



(Il999f ) at zs = 3 (panel [b]). 



The dot-dashed 
We find that 
models. 



ii ies in Figure 6 are the mean transmission models by iMadaul (119951 ). 



Madaul (119951 ) model almost always shows smaller transmissions than recent 



Madaul (119951 ) calculates the Lyman series regime from the equivalent width distri- 



bution of the LAF an d the Lyman cq i itinuu m regime from the column density distribution 



of the LAF and LLSs. 



Bershady et al. 



(1 19991) find that the Lyman series transmissions based 



on the equivalent width distrib ution of 



Madaul (119951 ) are smaller than those based on the 



column density distribution of iMadaul (1l995l ). and that the transr aissions i ncrea se further 



if the Doppler parameter b distribution is taken into accou nt, while IMadaul (119951 ) assumed 



a constant 6 = 35 km s ^. 



Tepper-Garcfa fc Fritzd (l2008l ) recently reproduces the Lyman 



a depression D\ for < z < 6 by their Monte Carlo simulation with the column density 



distribution of 



Madau ( 



alent with that of 



19951) not the equ ivalent width distribution, and their model is equiv- 



Bershady et al.l (119991) (but MC-N H model not MC-Kim model shown in 



Fig. 6). The smaller transmissions of IMadaul (119951 ) in the Lyman c ontinuum regi r ne for 
zs = I (panel [a]) are probably caused by the larger number of LLSs of ISargent et al.l (1l989l ) 
than ours (see Fig. 3). 



4.2 Lyman series and Lyman continuum transmissions 

Figure 7 shows transmissions averaged over the wavelength ranges blue- ward of Lyman a, /3, 
and 7 lines in the source rest-frame as a function of the source redshift. The case blue- ward of 
Lyman f3 line includes both of Lyman a absorption and Lyman (3 absorption. The case blue- 
ward of Lyman 7 line includes Lyman 7 absorption as well as Lyman a and (3 absorptions. 
Note that the transmissions shown in Figure 7 are wavelength-averaged ones, but those in 
Figure 6 are averaged over lines of sight. The filled circles with error-bars denote the median 
and the central 68% range in the distribution of the wavelength-averaged transmissions 
of 10,000 lines of sight for each source redshift. In the panels, we also show observational 
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Figure 7. Transmissions averaged over the wavelength ranges blue-ward of Lya (1070-1170 A; [a]), Ly/3 (980-1010 A; [b]), 
and Ly7 (955—970 A: [c]) in the source rest-frame. The horizontal axis is the source redshift. The filled circles with vertical 
error-bars are median and the central 68% range of the wavelength-averaged transmissions for 10,000 lines of sight generated 
in our Monte Carlo simulation. The open diamonds, circles, triangles, and asterisks are taken from Kirkman et al. (2007), 
Faucher-Giguere et al. (2008), Fan et al. (2006), and Songaila (2004), respectively. The open squares are the data obtained 
from Tepper-Garcia & Fritze (2008) but binned and corrected for the metal absorptions (the method by Tytler et al. 2004) by 
us. The cross mark is the estimation by us based on the observations by Steidel et al. (2001) (see appendix A). 



estimations for comparisons. We find an excellent agreement between our simulation and 
the observed data although a dispersion of the observed data is still found at ^ 3. This 
excellent agreement ensures the validity of our model for the Lyman series absorption. 

There are small upwards shifts of the filled circles (our model predictions) at the source 
redshift 2:3 = 0.2 relative to those at ^ 0.4 in the panels (b) and (c). These features 
are due to the lack of the Lyman a absorptions in the considering wavelength ranges. For 
example, to absorb the radiation of the Lyman f3 range (980-1010 A) of 2:3 = 0.2 by Lyman 
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Figure 8. Transmissions averaged over the wavelength range of the Lyman continuum (880—910 A) in the source rest-frame. 
The horizontal axis is the source redshift. The filled circles and vertical error-bars are median and the central 68% range of 
the wavelength-averaged transmissions for 10,000 lines of sight generated in our Monte Carlo simulation. The cross mark is the 
estimation by us based on the composite spectrum of 15 QSOs at {z) = 3.47 by Steidel et al. (2001) (see appendix A). 

a line, we would need absorbers with negative redshift. Since we consider the absorbers 
only at a positive redshift, the transmission in the Lyman f3 range shows such a small jump 
between = 0.2 and 0.4. The same is true for the Lyman 7 transmission. 

Figure 8 shows transmissions averaged over the wavelength range of the Lyman contin- 
uum in the source rest-frame as a function of the source redshift. This wavelength- averaged 
transmissions include all the Lyman series absorption and the Lyman continuum absorption. 
The filled circles and error-bars are the median and the central 68% range of the distribution 
of the wavelength-averaged transmissions of 10,000 lines of sight for each source redshift. We 
first find that the dispersion is quite larger than those of Lyman series wavelengths shown 
in Figure 7. The stochastic nature of the Lyman continuum absorption results in such a 
lar ge variance. We also show the observed estimation of the Lyman continuum transmission 



by 



Steidel et al 



(I2OOII ) in the figure. This is based on the composite spectrum of 15 QSOs 
at (z) = 3.47 (see appendix A for details). The value is compatible with our Monte Carlo 
simulation although it is somewhat smaller than the median of our simulation. In order 
to examine the validity of our model for the Lyman continuum absorption, we need other 
independent observational estimations of the Lyman continuum transmissions, which seems 
quite rare in the literature at the moment. We would encourage observational measurements 
of the Lyman continuum transmissions. The small jump in our model between the source 
redshift zs = 0.2 and 0.4 because of the same reason of Figure 7 (b) and (c). 
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Figure 9. Histograms of the logarithm of effective optical depths in the wavelength ranges blue-ward of Lya (1070—1170 A; 
[a]— [c]) and of the Lyman limit (880-910 A; [d]— [f]) in the source rest-frame. The source redshifts are noted in the panels. The 
thin solid histograms are the cases with only the LAF contribution (i.e. logjQ A^Hi/cm"^ < 17.2). The dotted histograms are the 
cases with the LAF and LLSs (log^o .'Vjji/cm"^ < 19.0). The dashed histograms are the cases without DLAs (logjQ A^jjl/cm"^ < 
20.3). The thick solid histograms are the cases with all types of absorbers. Each histogram consists of 10,000 lines of sight. 



4.3 Opacity distribution functions 

Figure 9 shows histograms of the logarithm of effective optical depths (Inrcff) corresponding 
to the wavelength-averaged transmission shown in Figures 7 and 8 (i.e. r^s = — ln(T)), of 
the Lyman a regime (the panels [a]-[c]) and of the Lyman continuum regime (the panels 
[d]-[f]) for three source redshifts. The thick solid histograms are the cases where all types of 
absorbers are taken into account. On the other hand, the dashed histograms are the cases 
without DLAs, the dotted histograms are the cases without DLAs and sub-DLAs, and the 
thin solid histograms are the cases of only the LAF contribution. 

The histograms of In TLyc, of all absorbers (thick solid ones in panels [a]-[c]) show a broad 
tail towards large opacity, whereas the histograms without DLAs do not show such prominent 
tails. Thus, we conclude that these tails are mainly caused by DLAs (A'^hi > 2 x 10^*^). Since 
the number density of DLAs increases along the redshift as shown in Figure 2, these tails 
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become more significant for liiglier redsliifts. Tlie = 5 liistogram of all absorbers seems 
largely modulated by the tail, whereas the fraction of lines of sight with InrLya > —1.0 for 
zs = 3 is 5%, and the fraction of lines of sight with InrLya > —2.5 for 2:3 = 1 is 1.4%. 
We also note that the effec t of LL Ss and sub-DLAs on TLya is small. This is already found 



by 



Tepper-Garcfa &: Fritzd (120081 ) with their Monte Carlo simulation and explained very 



well with the curve-of-growth theory. Since they did not include DLAs in their simulation, 
however, they could not find the tails found here. 

If we exclude DLAs, the histograms of In TLyo, (thin solid, dotted, and dashed histograms 
in panels [a]-[c]) seems Gaussian apparently. Note that the Gaussian distribution becomes 
parabola in the panels since the y-axis is shown in logarithmic scale. Suppose a null hy- 
pothesis that the histogram of In TLyo is Gaussian in order to examine the Gaussianity 
quantitatively. Based on the Kolmogorov-Smirnov (K-S) test, we cannot reject the null hy- 
pothesis with a significance level of 1% for the cases of only the LAF and of the LAF+LLSs 
at zs ^ 2. Thus, the distribution of TLy^ for these cases can be regarded as log-normal. 
For the = 1 cases, however, we can reject the null hypothesis with a significance level of 
0.01% by the K-S test. If we include sub-DLAs, Gaussianity of the distribution of In TLyo, is 
reduced, but few cases can be regarded as Gaussian (i.e. we cannot reject the null hypothesis 
with a significance level of 1% for the cases). 

We also find that the dispersion of InrLya decreases along the redshift. Note that the 
displayed range of the x-axis becomes narrower as the redshift becomes larger. Such a log- 
normal behaviour with a decreasing variance in In TLyo is found i n the observed probability 



distribution function of the transmitted flux of the QSOs' spectra (IBecker. Ranch, fc Sargent 



20071 ). On the other hand, the dispersion of transmissions shown in Figure 7 increases from 



z = to z = 3-4 and decreases towards higher redshifts (see also iTepper-Garcia fc Fritze 



20081 ). This difference between the transmission distribution and the logarithmic opacity dis- 



tribution is caused by logarithmic transformations from transmission to logarithmic opacity. 



Tepper-Garcfa fc Fritzd (120081 ) suggest that the distribution of the wavelength-averaged 
transmission for higher redshifts approaches Gaussian because of the Central Limit Theorem; 
sum of the contributions of a large number of absorbers. Indeed, the histograms of the 
wavelength-averaged transmission in the Lyman a regime of only the LAF case and of the 
LAF+LLSs case in our simulation could be regarded as Gaussian; we cannot reject the 
Gaussian hypothesis with a significance level of 1% for these cases at z^ ^ 3 based on the 
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K-S test. However, we also find with the K-S test that the distribution of InrLyo is closer to 
Gaussian than that of the wavelength-averaged Lyman a transmission. 

The histograms of Iutlc (Fig. 9 [d]-[f]) show a broad tail towards large opacity. Unlike 
luTLya, these tails are mainly produced by LLSs. Even in the LAF only case, we find such 
tails. Thus, the histograms of Iutlc are not Gaussian. These tails make the mode of the 
opacity always smaller than the mean. Thus, we have a large probability having a transmis- 
sion larger than that expected from the mean opacity model. We will come back to this point 
in section 5 where we discuss the detectability of the Lyman continuum of distant galaxies. 
We note that the cut-off at Iutlc ~ 4 found in the LAF+LLSs histograms of panels [e] and 
[f] corresponds to the optical depth for the upper limit of HI column density of LLSs (i.e. 
Nm ^ 1 X 10^^ cm~^ and tlc ^ 63 for a single LLS). 



4.4 Estimating Lyman continuum opacity from Lyman a opacity 

Can we estimate the Lyman continuum opacity from the Lyman a opacity? If it is possible, 
we may estimate the former opacity for an individual distant galaxy because we can estimate 
the latter opacity from the observed rest-frame 1000-1500 A spectrum individually. 

Figure 10 (a) shows the probability distribution in the plane of the logarithms of the 
effective Lyman a opacity (ri^ya) and the effective Lyman continuum opacity (tlc)- This 
is produced by 300,000 sets of (TLya,rLc): 30 source redshifts (0.2 - 6.0 w ith 0.2 interval) 



times 10,000 lines of sight. Contrary to the note by 



Shapley et al 



(120061 ) based on their 



Monte Carlo simulation of the intergalactic absorption, we find a good co r relatio n between 



InrLya and Iutlc; the correlation coefficient is 0.860. Since iShapley et al.l (120061 ) simulated 
the intergalactic absorption only for a single source redshift zs = 3.06, the dynamic range of 
TLya might be too small for them to find the correlation. Indeed, the probability density for 
Zs = 3 shown as the dotted contour in Figure 10 (a) is confined in a small range of XLyo, and 
elongated towards large tlc- Although LLSs mainly control tlc and the LAF does TLya, the 
LAF still has an effect on tlc- Thus, the correlation between tlc and TLya seems natural. 
However, we should confirm it observationally in future. Currently, we have one estimate 
cross mark i n the figure) based on the composite spectrum of 15 QSOs at {z) = 3.47 by 



Steidel et al 



(1200 ll ). which is consistent with our prediction. 
In Figure 10 (a), we also show the median and the central 68% range of the tlc distri- 
bution for 12 TLyQ, bins. Table 2 is the summary of them. Interestingly, we may estimate a 
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Figure 10. (a) Correlation between the effective optical depths in the wavelength ranges blue-ward of Lya (1070-1170 A; 
horizontal axis) and of the Lyman limit (880-910 A; vertical axis) in the source rest-frame. The contours show the probability 
density in this plane and the contour levels are 0.001, 0.003, 0.01, 0.03, and 0.1 from the outside; the solid contour is for 
2g = 0.2-6 and the dotted contour is for 23 = 3. The filled circles and the error-bars indicate the position of the median and 
the width of the central 68% range of the distribution in each bin of the horizontal axis (see Table 2). The cross mark is the 
estimation by us based on the QSO composite spectrum by Steidel et al. (2001) (see appendix A), (b) Difference of the medians 
and the central 68% ranges for two wavelength ranges of the Lyman continuum; filled symbols are the case of 880-910 A and 
open symbols are the case of 760-910 A. The open symbols are offset to the right by 0.1 dex for graphical clarity. 

statistically probable range of tlc and Lyman continuum transmission from TLya- We should 
note that the probable ranges of tlc depend on the wavelength range considered. In Fig- 
ure 10 (b), we show the difference of the medians and the 68% ranges for two cases of the 
wavelength range in the Lyman continuum: filled symbols are the case of 880-910 A and 
open symbols are the case of 760-910 A. For a wider wavelength case, we expect a larger 
optical depth (i.e. smaller transmission) because of a larger probability having a LLS in 
the considering wavelength range. Thus, we should care the wavelength range to estimate 
Tlc from Ti^ya- In Table 2, we have taken the wavelength range of 880-910 A in the source 
rest-frame. 



5 DETECTABILITY OF LYMAN CONTINUUM 



Although the observations of the Lyman continuum is very challengi n g bec ause of the 



intergalactic absorption as well as the interstellar one. 



Shapley et al. 



(120061) did detect 



the continuum from 2 out of 14 Ly man break galaxies at 2 ~ 3. On the other hand, 



Malkan. Webb, fc KonopackevI (l2008[ ) and 



Siana et al. 



(120071 ) did not detect the continuum 



from total 32 star-forming galaxies at z ~ 1. Combining these resul ts of direct observa tions 
with an indirect estimation from the ionising background intensity, llnoue et al.l (120061 ) sug- 
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Table 2. Median and central 68% width of the distribution of Imxc a-id corresponding Lyman continuum transmission as a 
function of In TLyc ■ 
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The wavelength ranges are 1070—1170 A for Lya and 880—910 A for the Lyman continuum in the source rest-frame. 

gested an evolution of the cosmic average escape fraction along the redshift: larger escape 
fraction at higher redshift on average. Such an evolution scena rio of the average escape frac- 



tion seems t o be supported 



2006 



2003). 



j y a rec ent cosmological simulation (IRazoumov fc Sommer-Larsen 



Shapley et al. 



(120061 ) showed that some galaxies do have a large escape frac- 
tion. If the number of the galaxies with a large escape fraction increases along the redshift, 
such galaxies may contribute to the cosmic reionisation significantly. Thus, understanding 
the properties of su ch galaxies is very in teresting. However, any common property in two 
detected galaxies of IShapley et al.l (120061 ) were not found. We may need much more number 
of galaxies showing a large escape of the Lyman continuum. Here, we discuss the detectabil- 
ity of such galaxies at z ^ 1 using the intergalactic transmission presented in the previous 
section. 

The observable flux density of the Lyman continuum is expressed as 



rpohs fp pobsT-iIGM 

-^LC ~ '^csc-Tuv-^LC ' 



(11) 



where T^esc is the escapin g flux density ratio of the Lyman continuum to the non-ionising 



ultraviolet introduced by 



Inoue et al 



(120061 ). is the observed non-ionising ultraviolet 



flux density, and Tlc'^ is the intergalactic transmission of the Lyman continuum. We may 
use the escape rati o T^esc as a proxy of the absolute escape fraction because they are the 
same order usually ( Inoue et al.l 2006 ). Estimating the absolute escape fraction is difficult in 
general because we need the dust attenuation correction and the intrinsic spec trum, while 



T^psr is just the observed flux density ratio corrected for the IGM opacity (see ilnoue et al. 



20061 for details). 

Suppose an observation in the Lyman continuum of distant galaxies with a limiting flux 
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Figure 11. Probability of having a line of sight with a transmission averaged over the Lyman continuum (880—910 A) in the 
source rest-frame larger than the value indicated in the horizontal axis. The lines from right to left are the cases of the source 
redshift from 1.0 to 5.8 with the redshift interval of 0.2. The cases of the source redshift of 1, 2, 3, 4, and 5 are indicated as 
the thick solid lines. 



density i^L™. In other words, we can detect a galaxy with its flux density in the Lyman 
continuum of Ff^ ^ with an enough significance. In this case, the IGM transmission 
towards the galaxy must satisfy 



IGM 

LC 



\ rplim 

^ -'lc = 



TT'lim 



fy T7obs 
'■^esc-f^UV 



(12) 



The quantity T^q means the minimum IGM transmission for a galaxy to be detected by an 
observation. Thus, the detectable probability of a galaxy is the probability that we have a 
transmission equal to or larger than Tl™. The detectable probability becomes larger for a 
fainter F^q (i.e. deeper observation), a larger T^esc (i-e. larger escape fraction), or a brighter 
(i.e. more luminous object). 

Figure 11 shows the probability of lines of sight with a Lyman continuum transmission 
Tlc'^ larger than the value indicated in the horizontal axis. This is just a cumulative prob- 
ability calculated from 10,000 transmissions produced by our Monte Carlo simulation for 
each source redshift ^s- We show 5 cases of zs = 1^5 as the thick solid lines. The case of 
2:3 = 6 is always zero probability in the range of the horizontal axis shown in the figure. 
We expect that the fraction of lines of sight not opaque against the Lyman continuum from 

= 3 (i.e. > 0.37 for tlc < 1) is 70%, and the fraction even for zs = 3.8 is 25%. 

Table 3 is a summary of the expected detectable probabil ity (%) with t h e assu mption 



of 7^, 



0.3 which is estimated from the two detections of 



Shapley et al 



fl2006h . When 



observing the Lyman continuum of 880-910 A in the source rest-frame, we find that 70- 
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Table 3. Detectable probability of a galaxy with a large escape fraction. 
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The detectable probability (%) of a galaxy with a large escape ratio of TZesc = 0.3 at the redshift = 1~5 when observing 
the wavelength range of 880-910 A in the source rest-frame. Seven cases of the limiting magnitude of the Lyman continuum 
observation are shown: mj^™'* = m'^y + 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, and 4.5, where m^'v observed magnitude in the 

non-ionising ultraviolet (~1500 A in the rest-frame) of the galaxy. We expect a larger probability for a deeper observation (i.e. 
larger "1^™'*) ^ more luminous object (i.e. smaller m'^y). Note that the detectable probability depends on the observing 
wavelength range (see text). 

80% galaxies with a large escape fraction of T^esc = 0.3 are detectable for zs = 3 and 
20-45% galaxies are detectable even for zg = A if the objects are 24.5 AB in UV and 
the limiting magnitu de of the Lyman continuum is 27.0-27.5 AB. In the observations of 



Shapley et al.l (120061 ). the limiting magnitude of the Lyman continuum and the average 



observed UV magnitude are 27.56 AB and 23.92 AB, respectively, and the average redshift 
is 3.06; th e detectable probabi lity of a galaxy with TZcsc = 0.3 in their sample is about 



90%. Since 



Shapley et al.l (120061 ) detected 2 out of 14 galaxies, therefore, we expect that the 



fraction of the galaxies with T^esc = 0.3 is about 16% in their sample. This suggests that 
there is a large variance of the escape fraction at z ~ 3, which is very interesting and whose 
origin should be understood in the context of the galaxy evolution in future. 

We note that the detectable probability in Table 3 may be the most optimistic case 
because the probability depends on the observing wavelength range. If we observe a wider 
wavelength range, the probability of having a LLS in the corresponding redshift range in- 
creases, then, the transmission and the detectable probability decreases. In addition, if we 
observe a wavelength range far from the Lyman limit (e.g., ~ 700 A), the probability affected 
by a LLS also increases, then, the detectable probability decreases (see also Figs. 4-6 and 
10). Therefore, we should take into account the observing wavelength range properly when 
estimating the detectable probability for future observations. 



6 CONCLUSIONS 

Motivated by recent attempts for determining the escape fraction of the Lyman continuum 
from distant galaxies, we have made a Monte Carlo simulation of the intergalactic absorption 
in order to model the Lyman continuum absorption properly. For this simulation, we have 
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derived an empirical distribution function of the intergalactic absorbers, Lyman a forest 

(LAF), Lyman limit systems (LLSs), and damped Lyman a systems (DLAs), from the 

observed statistics of the absorbers. In the distribution function of the absorbers, we have 

assumed a single functional form of the redshift distribution for all types of the absorbers, 

while the previous models of the intergalactic absorption assumed two different functional 

forms for the LAF and LLSs (and DLA). 

The transmission functions obtained from our simulation are very consistent with the 

I 1 

previous models except for IMadaul (119951 ) which predicts smaller transmissions than others. 
The Lyman series transmissions by our simulation excellently agree with the observational 
data, which ensures validity of our model. We have predicted the Lyman continuum trans- 
missions as a function of the source redshift. Although observational estimates of the Lyman 
continuum transmis sion are quite rare in the literature, one observational estimate based on 
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(I2OOII ) is consistent with our prediction. 
The distribution function of the Lyman a opacity for a source redshift seems log-normal 
with a tail towards large opacities. This tail is produced by DLAs. Since the number density 
of DLAs increases alon g the redshift, the signifi c ance of the tail in the distribution function 
increases. As found by iTepper-Garcia &: Fritzd (120081), we also find that the effect of LLSs 
on the Lyman a opacity is small. Along the source redshift, the mean of the log-normal part 
increases but the variance decreases. These features are found in recent observations. 

The distribution function of the Lyman continuum opacity shows a very broad tail to- 
wards a large opacity which is produced by LLSs. Unlike the Lyman a opacity, the effect 
of DLAs on the Lyman continuum opacity is small. Rarity of LLSs controlling the Lyman 
continuum opacity provides us with a chance to have a clean line of sight for z ~ 4; the 
probability of a clean (i.e. optical depth less than unity) line of sight at 900 A in the source 
rest-frame is about 70% for the source redshift z ~ 3 and about 20% for 2; ~ 4. 

A go od correlation betwee n the Lyman a opacity and the Lyman continuum one is found 
although IShapley et al.l (120061 ) noted no correlation in their simulation. This may be because 
a small dynamic range of the Lyman a opacity in their simulation. Based on the correlation 
we find, we may predict a statistically probable range of the Lyman continuum opacity as 
a function of the Lyman a opacity. This may be useful to estimate the Lyman continuum 
opacity from the Lyman a one for individual galaxy. 

Finally, we have predicted the detectable probability of a galaxy with a large escape 
fraction of the Lyman continuum based on our simulation. For example, the detectable 
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probability is 86% when the hmiting magnitude of the Lyman continuum observations is 
3.5 AB deeper than the non- ionising ultraviolet magnitude of the sample galaxy with an 
escape ratio T^esc = 0- 3 at z = 3. The small number fraction of the galaxies detected in the 
Lyman continuum by IShapley et al.l (120061 ) (2 out of 14 sample galaxies) suggests a large 
variance of the escape fraction a.t z = 3. This should be understood in the context of the 
galaxy evolution in future. 
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APPENDIX A: LYMAN a AND LYMAN CONTINUUM TRANSMISSIONS 
BASED ON QSO COMPOSITE SPECTRUM BY STEIDEL ET AL. (2001) 



We describe the method for estimating the transmissions in the 
continuum ranges based on the composit e spectrum of QSOs by 



yman a and the L yman 



Steidel et al 



(120011). The 



(120011 ). However, we assume a differ- 



method is essentially same as that used in lSteidel et al. 
ent intrinsic spectrum of QSOs based on the recent observations. The obtained transmissions 
are plotted in Pigures 7, 8, and 10. 

The wavelength ranges of the Lyman a and the Lyman continuum transmissions in 
Pigures 7, 8, and 10 are 1070-1170 A and 880-910 A, respectively, in the source rest-frame. 
We take 1100 A and 900 A as their representative values. In general, the IGM transmission 
at a wavelength x is defined as T^*^^ = F°^^/F^^^, where F°^^ and F™* are the observed and 
the intrinsic flux densities, respectively. The intrinsic flux density can be estimated from the 
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observed flux density at a wavelength free from the IGM absorption, that is, > 1216 A, if 
we assume an intrinsic spectrum. When we assume the reference wavelength of 1500 A, the 
Lyman a and the Lyman continuum transmissions are 

(Al) 



rplGM _ (-^110o/-^150o)obs 
(-FllOo/-P'l50o)mt 

and 

jnlGM _ {Fqoo/ Fi5oo)obs 
(-^900/ -^"1500)1111 

W e take the observed flux densities from the composite spectrum of QSOs by 



(A2) 
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( l200ll ). The sample consists of 15 QSOs at {z) = 3.47±0.14. The QSO compo site spectrum 
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i s con structed in order to compare the composite spectrum of LBGs from which 
(2001) detected the Lyman continuum of the sample LBGs. From Table 2 in 
fcoOlf ). we have F°}'^q = 1.86 ± 0.02 and F°}^^ = 0.50 ± 0.03 in the arbitrary unit of the flux 
density (per Hz). We obtain F°^qq = 3.4 from their value of (-Fi5oo/-^90o)obs = 6.8 ± 0.4. 

The intrinsic spectrum of QSOs is the largest source of the error in this estimation. The 
spectra of QSOs can be fit by a composite of power-laws. We here discuss the indices under 
the power-law assumption: F^, oc i/". Th e index in t h e ext reme-ultraviolet (including the 



Lyman continuum) is still controversial. 
-1.76 ± 0.12 from 184 QSOs at ^ ~ 1 (see also 



Telfer et all (2002) o 



btained a(500 - 1200 A) 



Scott et al. 



also 



Shang et al 



Zheng et al. 



19971 ). On the other hand, 



J2OO4 ) obtained a(630 - 1155 A) 



-0.56i[J:| from 85 AGNs at 2 ~ 0.2 (see 



2OO5I ). This con trast may be due to an evolution of the spectrum along 



the redshift and the luminosity (IScott et al.l 120041 ) , and also there may be the difficulty 
of the correction for the IGM absorption. We here assume two possible indices of a{< 
1100 A) = -1.8 ± 0.1 and -0.6 ± 0.3. For longer wavelength, Izheng et al.l Jl997n obtained 
a(1050 - 2200 A) = -0.99 ± 0.05. We here adopt a{> 110 A) = -1-0 ± 0.2 taking into 
account a large variance of the indices among QSOs (e.g., IShang et al.ll2005l ). Therefore, 
the intrinsic ratios of the flux densities adopted here are (-Fiioo/-Pi5oo)int = 0.73 ± 0.05, and 
(i^90o/i^i50o)int = 0.54 ± 0.04 or 0.69 ± 0.11. 

The observed r atios of the fl u x den sities measured in the composite spectrum of QSOs 
at (z) = 3.47 by ISteidel et al.l (I2OOII ) and the intrinsic ratios described above result in 
Tl^^ = 0.75 ± 0.05, and Tlf^ = 0.28 ± 0.03 or 0.22 ± 0.04. Finally, we summarise the 
Lyman continuum transmission as TIq^ = 0.25 ±0.06, where we put the sum of the internal 
and external errors as an error estimate. 



We should note one point here. In the procedure for obtaining the intrinsic spectrum 
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of QSOs, we need the IGM absorption correction. Then, we used the intrinsic spectrum 
to obtain T^'^^, Thus, it is hke a tautology although the intrinsic spectrum and T^^'^ are 
obtained from the different data and the different authors. However, we may consider that 
the obtained T^^'^ still have a meaning since we have assumed a wide range of the spectral 
indices for the intrinsic spectrum. For future measurements of T^G'^, we may use an inde- 
pendent intrinsic spectrum of QSOs, for example, based on an accretion disk model (e.g.. 



Kawaguchi. Shimura. fc Mineshigell200ll ). 



